Search for Cosmic Strings in CMB Anisotropics 
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ABSTRACT 

We have searched the lst-year WMAP W-Band CMB anisotropy map for evidence of cosmic 
strings. We have set a limit of 8 = 8ttGh/c 2 < 8.2 x 1CU 6 at 95% CL for statistical search for a 
significant number of strings in the map. We also have set a limit using the uniform distribution 
of strings model in the WMAP data with S = 8nGfi/c 2 < 7.34 x 1CT 5 at 95% CL. And the 
pattern search technique we developed here set a limit 6 = 8irGfj,/c 2 < 1.54 x 10 -5 at 95% CL. 



Subject headings: CMB anisotropy, cosmic strings 

1. Introduction 

Current theories of particle physics predict 
that topological defects would almost certainly 
be formed during the early evolution of the uni- 
verse [1]. Just as liquids turn to solids when the 
temperature drops, so the interactions between 
elementary particles run through distinct phases 
as the typical energy of those particles decreases 
with the expansion of the universe. When condi- 
tions favor the appearance of a new phase, the new 
phase crops up in many places at the same time, 
and when separate regions of the new phase run 
into each other, topological defects are the result. 
The detection of defects in the modern universe 
would provide precious information on events in 
the earliest moments after the Big Bang. Their 
absence, on the other hand, would force a major 
revision of current physics theories. The least of 
which would be to have the phase transitions be 
second or higher order. 

The potential role of cosmic topological defects 
in the evolution of our universe has interested the 
astrophysical community for many years. Com- 
bined theoretical and experimental work has led 
to the development of observational signatures in 



a variety of diverse data sets. These have, on one 
hand removed one major motivation of the pri- 
mary role in structure formation and on the other 
raised new motivations. These have also narrowed 
the range of allowable topological defects so that 
most viable models are likely to produce some va- 
riety of cosmic strings. 

CMB anisotropy power spectrum observations 
rule out topological defects as the primary source 
of structure in the universe [2,3]. Observations 
strongly favor adiabatic random fluctuations from 
something like inflation rather than the structure 
formation from topological defects (e.g. cosmic 
strings). These CMB results, while confirming in- 
flation, do not rule out lower level contributions 
from topological defects. Interest in cosmic strings 
has, in fact, been renewed with recent theoretical 
work on hybrid inflation, D-Brane inflation and 
SUSY GUTS. The idea that inflationary cosmol- 
ogy might lead to cosmic string production is not 
new; however, it has received new impetus from 
the brane world scenario suggested by superstring 
theory. A seemingly unavoidable outcome of brane 
inflation is the production of a network of cosmic 
strings [4], whose effects on cosmological observ- 



ables range from negligible to substantial, depend 
on the specific brane inflationary scenario [5]. A 
number of theory papers [12,13,14], anticipating 
that searches from cosmic strings will be negative 
and set significant limits, have begun developing 
modified theories so as not to produce them. At 
minimum these modified models must introduce a 
new field (also warping brane). 

Since a moving string would produce a step- 
like discontinuity in the CMB, it will cause tem- 
perature distribution deviate from Gaussian. We 
may be able to detect non-Gaussian aspects of 
temperature distribution provided sufficient res- 
olution. We search for strings in two ways: sta- 
tistical and pattern-discovery methods. Statisti- 
cal analysis determines how much the distribution 
of temperature fluctuation deviates from Gaus- 
sian distribution or what fraction of the fluctu- 
ations might be due to strings. In the second ap- 
proach, we can search for cosmic strings directly 
from temperature map by their distinctive pattern 
of anisotropy. The latter approach is much eas- 
ier and more straight forward when CMB signal 
is not contaminated seriously. These methods are 
distinctly different than a simple fitting to angular 
power spectrum [4]. 

2. Effect of Cosmic Strings on CMB 
2.1. Signal from a moving Cosmic String 

Consider a cosmic string with mass per unit 
length /i, velocity /3 and direction s both of which 
are perpendicular to the line of sight and is back- 
lighted by a uniform blackbody radiation back- 
ground of temperature T. Due to its angular de- 
fect 8 = 8ttG/i/c 2 , there is a Dopplcr effect of one 
side of the string relative to the other which causes 
a temperature step across the string [8,9] 



^ - SttG^/c 2 . 



(1) 



where 7 = 1/y 1 — (3 2 . This expression was gen- 
eralized for arbitrary angles between the string di- 
rection s, its velocity /3, and the line of sight n 

[6,7] 

AT 1 

— = SirG^Pn ■ x s)/c 2 . (2) 




Fig. 1. — Probability distributions for tempera- 
ture steps in equation (2). The solid curve rep- 
resents the probability with both string direction 
s and string velocity /3 having random directions 
while the flat dotted line assumes /3_Ls (sin cf> = 1) 
but the direction of /3 x s is random relative to the 
line of sight (cos 9 random) . 



2.1.1. Probability distributions for relevant pa- 
rameters 

ST/T in equation (2) is determined by three 
factors Gfi/c 2 , 7/3 and n ■ (/3 x s) which arise from 
different sources. Gfi/c 2 is related to symmetry 
breaking scale, we can assume 7/3 ~ 1 except near 
cusps where one can have 7 3> 1 and n ■ (/3 x 
s) depends on the combination of relative angles 
among /3, s and n. We can also assume that x s 
is random in direction to the line of sight n in 3-D 
space. If we denote /3 x s as u sin0, then 

n ■ (/3 x s) = h ■ u sin = cos 9 sin <fr . (3) 

The infinitesimal probability that cos 9 has an ar- 
bitrary value is 2-7T sin 9d9/4ir = ^dcos9 and thus 
its probability distribution is uniform in cos 9. If 
the string velocity /3 and its direction h are uncor- 
rected, the probability distribution for sin</> is y is 
j//yl — y 2 . Then, the probability distribution for 
k = n- (/3 x s) being an arbitrary value — 1 < k < 1 
becomes, substituting cos# = x, 
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Often string velocity and direction will be perpen- 
dicular to each other so that sin (f> ~ 1 and in this 
case the probability reduces to P(k) — \- These 
probability distributions are plotted in Figure 1. 

In a matter dominated universe the projected 
angular length of string in the redshift interval 
[z\, z?\ scales as ^fz{ — ^J~z~i [10]. The standard 
cosmological model (floE ~ 0.7, TO ~ 0.3) gives 
zi s ~ 1100 [10]. The apparent angular size of the 
horizon at the CMB last scattering redshift zi s is 
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(5) 

The average distance between strings is roughly 
dn/3. Thus the typical angular distance between 
the discontinuities on the sky to be of the or- 
der of Oh- The rough expected magnitude of 
the jumps in temperature are of order ST/T ~ 
13Gu/c 2 to 19G^/c 2 . With angular resolution of 
2' (2 arcmin) there could be sharp jumps up to 
about ST/T ~ 40G/x/c 2 . (i.e. peak range of tem- 
perature steps in a distribution of possible steps.) 
If the angular resolution is poorer, then the blur- 
ring effectively smooths the steps that the max- 
imum range of steps are at a somewhat smaller 
level. 

2.2. Statistical Fluctuations 

The signal in a microwave sky map will be made 
of several components: noise from the instrument, 
foreground signals (expected to be small away 
from the galactic plane and with strong sources 
punched out), CMB fluctuations from adiabatic 
(or appropriate) fluctuations, and potential signals 
from the strings. The signal observed at any pixel 
is then 

-Lpixel — 7V OlSGpixel t -t foreground ~r O-L CMB t Q± string • 

(6) 
Since the nature of signal is the superposition of 
random Gaussian signal (noise, STqmb) and non- 
Gaussian contribution (ST string ), there are two ba- 
sic forms of distribution functions (foreground will 
be removed from the beginning). 

Here and later on in this paper, we use two 
identical variables T and AT confusingly, i.e, T = 
AT = T abs oi ute - T (2.73K). 




Fig. 2. — Temperature Distribution of full 
sky from WMAP. Pixels with |AT| > lmK or 
|Galactic latitudc| < 10° are excluded because 
these signals are mainly due to the Galaxy or 
bright sources. The dashed curve is the best-fitted 
Gaussian with u — 23.QuK and a — 201.25[j,K. 
X 2 /DOF for the Gaussian fit is 4085.87/398 and 
the solid curve represents 2.6M-Gaussians fit with 
fi = 23.6fiK, ctq = 6682.22/ii4T and ocmb = 



80.25/jK. 
435.76/398. 



X /DOF for 2.6M-Gaussians fit is 



(l)Gaussian signal 

Major portion of CMB signal obeys Gaussian dis- 
tribution, 



h(T) 
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(7) 



2 2 / i 2 

a = a /n % + a CMB 



where o§ and ocmb are variances from noise and 
CMB fluctuations each, rii is the number of obser- 
vations for i t h pixel. We make use of the informa- 
tion that the noise is Gaussian (actually Gaussian 
per observation and thus variance inversely pro- 
portional to the number of observations per pixel) 
and the intrinsic CMB fluctuations are gaussianly 
distributed. 

(2) Non- Gaussian signal from strings 
When a straight moving string is added to a re- 
gion, we can expect a temperature distribution 
that is the sum of two Gaussians rather than sin- 
gle Gaussian because of the blue- and redshift by a 
transversely moving string. Then the probability 
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Limit G/i/c 2 GUT symmetry brcakinj 
xlO 6 scale r, SB (xl0 16 GcV) 



Signal 



Total 


201.25 


7.37 xlO" 5 


2.93 


Total - Noise 6.7mK/obs 


80.26 


2.94xl0~ 5 


1.17 


Total - Noise - Adiabatic CMB 


22.49 


0.82 xl0~ 5 


0.33 


Uniform Distribution of Strings 


57.82 


7.34xl0~ 5 


2.92 


Pattern Search 


N/A 


1.54xl0- 5 


0.61 



2.05 
1.30 
0.69 
2.05 
0.94 



Table 1: Variance Limits on Cosmic Strings. All the values presented are at 95% CL. The model 'Uniform 
Distribution of Strings' on the fourth row is the distribution given in the equation (18). The GUT 
symmetry breaking scales, tisb, in the last column are calculated using the relation tjsb ~ m p ; -v/G/i/c 2 
with m p i = 1.2 x 10 19 GeV. 



distribution in equation (7) is modified to 
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where p — Nu U e I 'N to tai , the ratio blucshiftcd pix- 
els to total pixels, q = 1 — p and ST/ 2 is half of the 
effective height of step across the string as given 
in equation (2). If we denote v and a b s as mean 
and standard deviation calculated from observa- 
tion, they are related to jjl and a as follows 
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= a 2 +4p(l-p)(ST/2) 2 . 


(10) 



Here we used an approximation on the integra- 
tion interval as — ImK < T < ImK — ► — oo < 
T < oo. It is a good approximation because, a ~ 
200fiK and thus ImK ~ 5cr, the variance integral 
in equation (10) over [— 5cr, ha] covers 99.9985% of 
that over [-co, oo]. This approximation is also im- 
plied in equation (19). Temperature distribution 
of sky may be explained with appropriate combi- 
nations of fi(T) and f2(T). All the models intro- 
duced in the following sections are based on these 
two primary patterns of temperature distribution. 

2.2.1. Variance: Quadratic Estimator Test for 
String Contribution 

Assuming that the signals are all statistically 
independent, we can find the variance in the map 
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foreground 



} CMB 



string 



(ii) 



We remove the significant foreground contribution 
by dropping pixels with | Galactic latitude| < 10° 
or |AT| > ImK. The first condition is to ex- 
clude galactic area where non-CMB signal is dom- 
inant. We lose 217 more pixels by imposing the 
condition |AT| < ImK, 4 pixels are less than 
— ImK and 213 pixels are greater than ImK. 
Since Gaussian tail probability allows ~ 1.7 pixels 
in | AT | > ImK region and they mostly form clus- 
ters, we can drop those pixels as they are from un- 
usual bright sources. We introduce a temperature 
distribution based on the Gaussianity of signal for 
each pixel in WMAP lst-year data and calculate 
o'cmb which contains all the contributions except 
instrumental noise, 



f2.6M{T) — 



1 E 



1 



'CMB 



(12) 



where N is the total number of pixels, 2,598,695, 
i is pixel number with maximum resolution of 
WMAP lst-year data, a 1 is variance due to noise 
per measuremet and n» is effective number of mea- 
surements for %th pixel. The distribution /2.6m(T) 
in equation (12) has minimum x 2 for o~cmb = 
80.25t?;65A^ and a = 6682.22l| 1 ^ 8 / u^ and thus 
a s = 200.68±g;|^iir each with 95% CL. This ob- 
servation in turn allows an upper limit of order 
< 

^string ~ 



8irGn/c 2 - 200^^/2.73^ 



200/iK and thus providing a limit of 
73.5 x 10~ 6 . 
We obtain a better upper limit, if we know the 
mean contribution of instrumental noise and other 
signals. According to the WMAP lst-year data re- 
lease, the sum over pixels of (number of measure- 
ments per pixel) -1 for 2,598,695 pixels included 
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Fig. 3. — Temperature distribution of full sky ex- 
cluding the pixels with |Galactic latitude| < 10° 
or | AT | > lmK. Solid curve represents the best- 
fitted curve of uniformly distributed strings model 
and the dotted curve shows actual WMAP 1st- 
year data. 

here is 1968.96. Thus the effective variance due to 
noise is 

oLise = K 2 £-J-W 

i 
= ( 6682 - 22 -3-) 2 ^^ 

= (183.93i^ 8 A^) 2 , 95% CL (13) 

where rn is the effective number of measurements 
for i th pixel. 

Under the assumption that the signals other 
than noise and from strings are negligible, 
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'string — a S 



.25 



+0.05 
-1.65 



ixKf 



95% CL. 

(14) 

And we are able, by taking the instrument noise 
into account, to set a lower upper limit. 

Taking into account the CMB anisotropy spec- 
trum is characteristic of anisotropies that arise 
from stochastic, Gaussianly-distributed, adiabatic 
primordial perturbations. We decompose the vari- 
ance (Tq MB computed in equation (14) into two 
contributions, one from adiabatic fluctuation and 
one from strings, a 2 CMB = o 2 adlabatic + a 2 trlng . 
Since cf cmb = ^(2Z + l)C;.B; 2 /47r and string con- 
tribution in C\ is ~ 10% [4] where Bi is beam 
filter function, we estimate the upper limit of 



Fig. 4. — x distribution for the uniformly dis- 
tributed strings model. Xmin ~ 435.43 with 398 
degrees of freedom. Parameters span < ST < 
200.3 iiK and 192.1 m# < a < 200.8 fj,K at 95% 
CL. Here a 1 represents the total variance except 

2 

string contribution, a 1 = ^ J^ 1/nj + o'cmb- 
The outer contour is for 95% CL which is almost 
overlapped with the 68% CL contour in vertical 
direction. 



String ~ ^CMb/^ = (25.381?$ ^K) 2 , thus 

< a strmg < 25.38t?:oV^ 95% CL. (15) 

We can also calculate cr 2 diabatic directly from the 
relation between cr CMB an d the power spectrum 
coefficients Ci . Using the Ci 's of best-fitted mod- 
eled cosmological parameters (ft\ = 0.73, fib = 
0.046, Q cdm = 0.224), B ; 's from WMAP lst-year 
data and convolving with finite pixel size, we have 

' ^(2l + l)QB? 



2 
"CMB 



47T *-^ 

I 



= (79.085 ±0.376^) 2 , 95% CL.(16) 

for beam and finite pixel size. Assuming that 
a CMB m equation (16) is solely from adiabatic 
fluctuation, we have even lower upper limit on 

^string ; 



string 



< ai 



adiabatic 



= (201.25/iiT) 2 - (183.85m-?0 2 
= (22A9fj,K) 2 , 95% CL. 



(78.709^.?0 2 

(17) 



Both of these last two are model dependant as one 
could refit all cosmological parameters to include a 
small addition of string. Presumably that is what 
was done in [4]. 
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Fig. 5. — Generating the vector field L. Pixels left demoted to middle pixels, middle 
right: defining gradient of temperature VT and L which lies along temperature step. 



2.2.2. Statistical Test for Temperature Step Ex- 
pected from Random Strings 

For the anticipated distribution of strings, one 
would expect a random distribution of tempera- 
ture steps which has roughly equal probability be- 
tween the plus and minus the maximum temper- 
ature step amplitude. We also fitted to a distri- 
bution that was a Gaussian for the other signals 
convolved with a uniform distribution of temper- 
ature steps. Let 8Tq = 87rG/i7/?/c 2 To be the char- 
acteristic value for a string, then the actual effect 
of temperature step left on the measured CMB is 
kSTo/2 where — 1 < k = cos 9 smcf> < 1 from the 
equation (4) and Tq = 2.73K. Assuming sin</> ~ 1 
i.e direction of a string and its velocity are most 
likely perpendicular, we have P(k) ~ 5- Then, we 
obtain the distribution function of temperature 



m^h^E//"' 



(T,k5T /2,a)dk (18) 



where f 2 (T, k5To/2, a) is the two-Gaussians form 
with p = 0.5 and 6T = k5To of i-th Gaussian 
in equation (12). Then, the apparent variance 
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uniform 



for this distribution becomes, in terms 



of a and STq , 
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(T - nff T (T)dT 
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r^ + ^f (19) 



where the first 2 terms are directly from the equa- 
tion (12) and the last term is the contribution 



from strings. From the x 2 distribution in Fig- 
ure 4, we have 192.1 ^.K < a < 200.8 fj,K and 
< ST < 200.3 fj,K at 95% CL. Thus, tempera- 
ture variation by a moving string can be limited 
to, in terms of deficit angle, 



<^£ = ^<7.34xl0- 5 ,95%CL. 



2.3. Search for Temperature Steps 



(20) 



One can directly search for temperature steps 
from the CMB sky map by their topological con- 
figuration or pattern. We define an algorithm to 
search for CMB temperature steps produced by 
cosmic strings in a background of CMB adiabatic 
fluctuations and (at present even more dominant) 
receiver noise. (For the WMAP lst-year data re- 
lease the signal-to-noise ratio for an average W- 
band pixel is about 0.5.) As a smoothing pro- 
cess, we demote 2 2 ™ pixels of maximum resolu- 
tion to one pixel with representative temperature 
value (average) assigned. This demoted pixel can 
be obtained by repeating demotion process shown 
in Figure 5 until we reach desired value of sig- 
nal/noise. Since as ~ 201/iZf and a n oise ~ lMfiK 
for the WMAP lst-year data, it is reasonable to 
take 16 pixels demotion or more because we have 



r (16) 



a x s "' - 104fj,K and a noi ' se = a noise /VlQ ~ AQ^K, 



thus the contribution from noise becomes subdom- 
inant. Then a vector can be derived out of 4 neigh- 
boring demoted pixels. VT is defined as 



(VT), 
(VT), 



(T 2 +T 4 -T 1 -T3) (21) 
(T 3 + T i -T 1 - T 2 ) . (22) 
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Fig. 6. — L field for the full sky with 64-pixel demotion. The center of the map is (0 = 0°, <f> = 0°) 
in Galactic coordinates. Each arrow represents local temperature step and the temperature steps 
greater than 100/j,K are set to 100/iK. The pattern along the equator is due to the Galaxy. 



We build up a vector field for the full sky which 
contains 12 x 2 w ~ 2n VT's where 2 2n is the de- 
motion factor. For each VT, the perpendicular 
direction to VT represents the direction of step 
discontinuity. So if there is a consistent elongated 
line of discontinuity we may interpret it as the 
signal due to moving string. After covering full 
sky with VT-field, we rotate the vector field —ir/2 
to define L-field so that each arrow lies along the 
local isothermal line with lefthand side of arrow 
being higher temperature. Figure 6 shows the L- 
field map for the entire sky. The length of each 
arrow represents the magnitude of temperature 
gradient on that point and any temperature step 
greater than 100/j.K is set to 100/iK to make steps 
visible which are less than 100/xif . It shows a long 
coherent structure along the galactic plane, this 
is because there is a steep temperature rise ap- 
proaching the galactic plane. 



2.4. Evaluating String Pattern 

We define the connectedness of two neighbor- 
ing temperature steps with two smoothness con- 
ditions for the heights of steps and the curve that 
links neighboring L's. 

(1) Component Definition of Connectedness 
We assume that the temperature distribution of 
WMAP lst-year data is approximately Gaussian 

(Figure 2), f T (T) = (ctsv^)- 1 ^^^ 2 with 
ji ~ 25fiK and as — 200/iK. Starting with this 
temperature distribution fr{T), after taking n- 
pixcl demotion, we can derive the distributions of 
Lj™ anci L 1( , 



J compy^ x ) 



1 



2" 



(23) 



where a 2 is the variance for temperature distri- 
bution for n-pixel demoted pixels for example, 
CTi6 — 104/iK. L y also has almost same distri- 
bution as equation (23), so we use equation (23) 
for both L x and L y . This gives the probability 




Fig. 7. — Relative angles defined for connected- 
ness. £>i is a unit vector along the line that con- 
nects Lj_i and Li. di.i+\ is defined by the angle 
between Di and £>j+i, Oi,i+i = cos (Di • -Di+i) 
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distribution function for change of x-component, 

L± x — L, L, , 



/a.(A x ) = 



1 



2a n ^ 



(24) 



and the same function for y-componcnt. Since a 
pattern formed by moving string should have rel- 
atively constant height of step along the curve of 
pattern, we can impose a condition for connected- 
ness as 

\A X \ < AL C , \A y \ <AL C (25) 

then, the probability that two adjacent vectors Li 
and L i+ i meet these conditions is given by 



Pi 



»,»+! 
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AL C 

erf(-= — ) 



f Aas (A x )dA x 



(26) 



where AL C is the maximum value allowed for A a 
and A y . 

We set another condition for connectedness 







i,i+l 



COS \Di • A+l) < 



(27) 



for a sequence to avoid too sharp turns (Figure 7) 
where C is the maximum angle allowed for Qi.i+\ 
to claim Li and Li+i are connected. Figure 8 
and 9 show examples of patterns which comply 
the definition described in equations (25) and (27) 
found in the WMAP lst-year data. Here, we set 
AL C = a w = 104fj,K and 9 C = tt/6. 
(2) Likelihood of sequence 



Fig. 8. — A long temperature discontinuity pat- 
tern found at (I = 45°, b = 43°) in Galactic coor- 
dinates. The vector field is derived from 16-pixel 
demoted temperature sky map. This sequence has 
its maximum likelihood defined in the equation 
(28) at ST = 12/j.K. 



Given that a sequence of temperature steps de- 
fined in the previous paragraph, we can estimate 
its likelihood for a signal due to a moving string. If 
a long temperature step is formed solely by a mov- 
ing string, the L's assigned on the curve should be 
tangent to the curve. After allowing contamina- 
tion by noise, adiabatic fluctuation or by other 
possible sources, each L on the curve will be off 
from local tangent. But if the contamination is not 
overwhelming, there should be still a bias seeded 
by string. We define the bias of a sequence with 
iV-connected arrows quantitatively in terms of rel- 
ative likelihood function as follows 



A' 



42(*n = He-(£M)/^ 



(28) 



i=i 



right — handed curve 

straight line 
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where STi is a tangent vector at ith grid point on 
the curve at which Li is located and \6Ti\ = ST. 
<t>i is a phase factor defined at « t h grid point on 
the curve to give correct direction of string veloc- 
ity on the curve. If the sequence turns right (left) 
locally at i t h grid, then (pi = +1(— 1) and it is zero 
when the sequence is straight at that point. On 
both the head and tail of a sequence, the phase 
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Fig. 9. — A short pattern with 16-pixel demotion 
from WMAP lst-year data located at I - 127°, 
b = 20° in Galactic coordinates. ST opt for this 
pattern is 0/iK, i.e, no preferred bias is found 
which is against its apparent feature. This is be- 
cause the temperature step vectors on the left-half 
of the curve have wrong direction relative to ex- 
pected string velocity at this part while those on 
the right-half arc aligned correctly and these two 
opposite contributions cancel out the bias. 



factors are set to zero. This is an approximate 
prescription to describe realistic model of string 
motion. When Lj's are perfectly coincident with 
STi, then £ bias ($T) become 1 and if Li's are off by 
either direction or magnitude or both from STi's, 
then £ bias ($T) decays exponentially. So, if there 

is nonzero ST that gives maximum C bias (ST) , then 
it is relatively more likely that there is a constant 
temperature step with height ST imbedded in the 
sequence. Figure 10 shows the comparison of re- 
sults between actual data (WMAP lst-year data) 
and a simulated white noise. We can estimate 
from the curve the height of temperature discon- 
tinuity in equation (1), 



ST<42fiK, 95%CL 



(29) 



which is equivalent to symmetry breaking scale 
tjsb ~ 0.94 x 10 16 GeV. 

3. Conclusion 

We have investigated WMAP lst-year data to 
search directly or set limit on cosmic strings. For 



Fig. 10. — Distribution of ST op t at which likeli- 
hood function defined in equation (28) becomes 
maximum for a sequence. Each data set takes 
16-pixel demotion. We included sequences which 
contain five or more arrows connected. There are 
many more sequences in the actual data than in 
the simulation which has identical size, mean and 
variance but white noise. 



statistics of full sky map, we set the limit of contri- 
bution to variance by the strings as < a string < 
22. 5 /iK which gives upper limit on deficit angle 
< S = 8ttG^/c 2 < 0.82 x 10~ 5 . And we set 
up a model with a random distribution of strings 
with random orientation with which the relation 
between a string and height of temperature step 
8ttG^,j(3/c 2 can be found. A limit on deficit angle 
is obtained from the model of uniform distribution 
of strings, < 6 = 8irGn/c 2 < 7.34 x 10~ 5 . This 
corresponds to symmetry breaking energy scale 

^-2x 10 16 GeV. 



■qsB ~ m p i * 

We developed a pattern search algorithm that 
can visualize the landscape of CMB temperature 
variation of the sky. There were some fairly long 
temperature rows but we didn't find any com- 
pelling pattern of cosmic strings predicted by the- 
ory. Instead, by considering the distribution of 
heights of temperature discontinuities, we roughly 
estimated < S = 8irG(i/c 2 < 1.54 x 10~ 5 or 
equivalently iqsB ~ 0.94 x 10 16 GeV. 

The precision of data is yet to be refined and we 
expect WMAP 2nd-year data will provide much 
better chance to pin down the effects of topolog- 
ical defects on cosmic microwave background ra- 



diation and at that point refined analysis will be 
appropriate. 
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